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Abstract 

We consider the complexity of Delaunay triangulations of 
sets of points in IR^ under certain practical geometric 
constraints. The spread of a set of points is the ratio between 
the longest and shortest pairwise distances. We show that 
in the worst case, the Delaunay triangulation of n points 
in IR^ with spread A has complexity 0(min{A^ , nA, n^}) 
and OlminlA'^.n^}). For the case A — 0{^/n), our lower 
bound construction consists of a uniform sample of a smooth 
convex surface with bounded curvature. We also construct a 
family of smooth connected surfaces such that the Delaunay 
triangulation of any good point sample has near-quadratic 
complexity. 



1 Introduction 

Delaunay triangulations and Voronoi diagrams are used 
as a fundamental tool in several geometric application 
areas, including finite-element mesh generation [16, 21, 
27, 28], deformable surface modeling [15], and surface 
reconstruction [1, 3, 4, 5, 12, 26]. Many algorithms in 
these application domains begin by constructing the 
Delaunay triangulation of a set of n points in IR^. 
Delaunay triangulations can have complexity Ofn^) in 
the worst case, and as a result, all these algorithms 
have worst-case running time D.i'n?]. However, this 
behavior is almost never observed in practice except 
for highly-contrived inputs. For all practical purposes, 
three-dimensional Delaunay triangulations appear to 
have linear complexity. 

One way to explain this frustrating discrepancy be- 
tween theoretical and practical behavior would be to 
identify geometric constraints that are satisfied by 
real-world input and analyze Delaunay triangulations 
under those constraints. These constraints would be 
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similar to the realistic input models such as fatness or 
simple cover complexity, which many authors have used 
to develop geometric algorithms with good practical 
performance [8, 29]. Unlike these works, however, our 
(immediate) goal is not to develop new algorithms, but 
rather to formally explain the good practical perfor- 
mance of existing code. 

Dwyer [20, 19] showed that if a set of points is 
generated uniformly at random from the unit ball, its 
Delaunay triangulation has linear expected complexity. 
Golin and Na [24] recently derived a similar result for 
random points on the surface of a three-dimensional 
convex polytope. Although these results are encourag- 
ing, they are unsatisfying as an explanation of practical 
behavior. Real-world surface data generated by laser 
range finders, digital cameras, tomographic scanners, 
and similar input devices is often highly structured. 

This paper considers the complexity of Delaunay 
triangulations under two types of practical geometric 
constraints. First, in Section 2, we consider the worst- 
case Delaunay complexity as a function of both the 
number of points and the spread — the ratio between 
its diameter and the distance between its closest pair. 
For any n and A, we construct a set of n points with 
spread A whose Delaunay triangulation has complexity 
a{min{A3,nA,n2}). When A = 8(Vn), our lower 
bound construction consists of a grid-like sample of a 
right circular cylinder with constant height and radius. 
We also show that the worst-case complexity of a 
Delaunay triangulation is OfminjA'^.n^}). 

An important application of Delaunay triangulations 
that has received a lot of attention recently is surface 
reconstruction — given a set of points from a smooth 
surface L, reconstruct an approximation of Z. Several 
algorithms provably reconstruct surfaces if the input 
points satisfy certain sampling conditions [4, 5, 12, 26]. 
In Section 3, we consider the complexity of Delaunay 
triangulations of good samples of smooth surfaces. 
Not surprisingly, oversampling almost any surface can 
produce a point set whose Delaunay triangulation has 



quadratic complexity. We show that even surface data 
with no oversampling can have quadratic Delaunay 
triangulations and that there are smooth surfaces where 
every good sample has near-quadratic Delaunay com- 
plexity. We also derive similar results for randomly 
distributed points on non-convex smooth surfaces. 

We will analyze the complexity of three-dimensional 
Delaunay triangulations by counting the number of 
edges. Two points are joined by an edge in the Delaunay 
triangulation of a set S if and only if they lie on a sphere 
with no points of S in its interior. Since every vertex 
figure is a planar graph, Euler's formula implies that a 
Delaunay triangulation with n vertices and e edges has 
at most 2e — 2n triangles and e — n tetrahedra. 

In the interest of saving space, several straightforward 
but tedious calculations are omitted from this extended 
abstract. 

2 Sublinear Spread 

We define the spread A of a set of points (also called the 
distance ratio [17]) as the ratio between the longest and 
shortest pairwise distances. In this section, we derive 
upper and lower bounds on the worst-case complexity 
of the Delaunay triangulation of a point set in IR , as a 
function of both the number of points and the spread. 

If the spread takes its minimum value 0[n^/^], the 
points are packed into a tight lattice, and the Delaunay 
triangulation has only linear complexity. On the other 
hand, all known examples of point sets with quadratic- 
complexity Delaunay triangulations have spread Ofn). 
Thus, it is natural to ask how the worst-case complexity 
of the Delaunay triangulation changes as the spread 
varies between these two extremes. The spread of a set 
of points is loosely related to its dimensionality. If a set 
uniformly covers a bounded region of space, a surface 
of bounded curvature, or a curve of bounded curvature, 
its spread is respectively 0(n^/^), 0(n^/^), or 0(n). 
The case of surface data is particularly interesting in 
light of numerous algorithms that reconstruct surfaces 
using a subcomplex of the Delaunay triangulation. We 
will discuss surface reconstruction in more detail in the 
next section. 

2.1 Lower Bounds 

The crucial special case of our lower bound construction 
is A = Q{y/n]. For any positive integer x, let [x] 
denote the set {1,2, ...,x}. Our construction consists 
of n evenly spaced points on a helical space curve: 

t . t 



See Figure 1. As we show below, the Delaunay 
triangulation of has complexity n(n^/^). Note 
that S is a grid-like uniform £-sample of a right circular 
cylinder, where £ — Q(^y^ /n). By adding additional 
points on two hemispherical caps at the ends of the 
cylinder, we can extend S into a uniform sample of 
a smooth convex surface with bounded curvature and 
constant local feature size. 




Figure 1. A set of n points whose Delaunay triangulation has 
complexity Cl{n^''^) 



Let h.cc(t) denote the helix (at, cost, sin t), where 
a > is a fixed parameter called the pitch. Using ele- 
mentary trigonometric identities and matrix operations, 
we can simplify the insphere determinant for five points 
on this helix as follows. 
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We obtain the surprising observation that changing the 
pitch a of the helix does not change the combinatorial 
structure of the Delaunay triangulation of any set of 
points on the helix. (More generally, scaling any set 
of points on any circular cylinder along the cylinder's 
axis leaves the Delaunay triangulation invariant.) Thus, 
for purposes of analysis, it suffices to consider the case 
a= 1. Let h,(t) = hi (t) = (t, cost, sint). 

Our first important observation is that any set of 
points on a single turn of any helix has a neighborly 
Delaunay triangulation, meaning that every pair of 
points is connected by a Delaunay edge. For any real 
value t, we define the bitangent sphere |3(t) as the 
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unique sphere passing through h.(t] and h.(— t] and 
tangent to the hehx at those two points. 

Lemma 2.1. For any < t < n, the sphere |3(t] 
intersects the helix h, onJj at its two points of tangency. 

Proof: Symmetry considerations imply that the bitan- 
gent sphere must be centered on the y-axis, so it can be 
described by the equation x^ + (y—a)^+z^ — for some 
constants a and r. Let y denote the intersection curve 



Theorem 2.3. For any n, there is a set of n points 
in IR^ with spread \/n whose Delaunay triangulation 
has complexity n(n^/^). Moreover, this point set is 
a uniform sample of a smooth convex surface with 
constant local feature size. 

We can generalize our helix construction to other 
values of the spread A as follows. 

Theorem 2.4. For any n and A = Q(n'/^), there is a 



of (3(t) and the cyhnder y^+z^^^. Every intersection °^ ^ P"'""^' •'^ ^ ^^'^^ ^P^^^^ ^ Delaunay 



point between (3(t) and the helix must lie on y. If we 
project the helix and the intersection curve to the xy- 
plane, we obtain the sinusoid y = cosx and a portion of 
the parabola y = y(x) = (x^ — + + 1 )/2a. These 
two curves meet tangentially at the points (t, cost) and 
(— t, cost). 




Figure 2. The intersection curve of the cylinder and a bitangent 
sphere projects to a parabola on the xy-plane. 

The mean value theorem implies that y(x) — cosx at 
most four times in the range — tt < x < tt. (Otherwise, 
the curves y" — — cosx and y" — t"(x) — 1/a would 
intersect more than twice in that range.) Since the 
curves meet with even multiplicity at two points, those 
are the only intersection points in the range — tt < x < 
TT. Since y{x) is concave, we have y(±tt) < cos±7T = 
— 1 , so there are no intersections with |x| > tt. Thus, the 
curves meet only at their two points of tangency. □ 

Corollary 2.2. Any set S of n points on the helix h,(t) 
in the range — tt < t < tt has a neighborly Delaunay 
triangulation. 

Proof: Let p and q be arbitrary points in S, and let |3 
be the unique ball tangent to the helix at p and q. By 
Lemma 2.1, |3 does not otherwise intersect the helix 
and therefore contains no point in S. Thus, p and q are 
neighbors in the Delaunay triangulation of S. □ 

We can now easily complete the analysis of our 
helical point set S^. Lemma 2.1 implies that every 
point in is connected by a Delaunay edge to 

every other point less than a full turn around the 
helix h^yYj:^(t), and each full turn of the helix contains 
[2tty^J points. Thus, the number of edges in the 
Delaunay triangulation of is at least Ztttt^/^ — 0(n). 



triangulation has complexity n(min{A^ , nA, n^}) . 

Proof: There are three cases to consider, depending on 
whether the spread is at least n, between ^yr\. and n, 
or at most -^n. The first case is trivial. For the case 
-^/n < A < n, we take a set of evenly spaced points on 
a helix with pitch A/n: 
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Every point in Sa is connected by a Delaunay edge to 
every other point less than a full turn away on the helix, 
and each turn of the helix contains Ci[A) points, so 
the total complexity of the Delaunay triangulation is 

a(nA). 

The final case n'^^ < A < y/ri is slightly more 
complicated. Our point set consists of several copies 
of our helix construction, with the helices positioned at 
the points of a square lattice, so the entire construction 
loosely resembles a mattress. Specifically, Sa is the set 
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where r and w are parameters to be determined shortly. 
This set contains n — w^r points. The diameter of Sa 
is 0[w) and the closest pair distance is 0[]/\/r), so 
its spread is A = @(wy/r). Thus, given n and A, we 
have w = Q(n/A^] and r = Q(A^/n^). Straightforward 
calculations imply that for all t < tt/4 and a < 1, the 
bitangent sphere (3a(t) has radius less than 2. Since 
adjacent helices are separated by distance 2, every 
point in Sa is connected in the Delunay triangulation 
to every point at most half a turn away in the same 
helix. Each turn of each helix contains D.(^/r) points, 
so the Delaunay triangulation of Sa has complexity 
a(nVr) =a(A3). □ 

2.2 Upper Bounds 

Let B be a ball of radius R in IR^ , and let b i , , ba , . . . 
be balls of radius at least r, where 1 < r < R. Our upper 
bound proof uses the following geometric properties of 
the 'Swiss cheese' C = B \ IJ^ bi. See Figure 3. 
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Figure 3. Swiss cheese (in IR^) 

Lemma 2.5. Tie surface area of C is 0(R^/r]. 

Proof: The outer surface 9C n 9B clearly has area 
O(R^) — 0(R^/r), so it suffices to bound the surface 
area of the 'holes'. For each i, let Ht = B n 9bi be the 
boundary of the ith hole, and let H = IJ^ = 9C \ 9B. 
For any point x g H, let Sx denote the open line segment 
of length r extending from x towards the center of 
the ball bt with x on its boundary. (If x lies on the 
surface of more than one bt, choose one arbitrarily.) 
Let S — UxgH union of all such segments, and 

for each i, let Si — Uxsh Each St is a fragment of 
a spherical shell of thickness r inside the ball b^. See 
Figure 4. 




Figure 4. Shell fragments used to bound the surface area of C. 

For each i, we have (after some tedious calculations) 

vol(Si) = - ^—^^2^^ area{Hi] > ^ area(Hi), 

where ti > r is the radius of hi. The triangle inequality 
implies that Sx and Sy are disjoint for any two points 
x,y G H, so the shell fragments Si are pairwise disjoint. 
Finally, since S fits inside a ball of radius R + r < 2R, 
its volume is O(R^). Thus, area(H) — ^^^area(Hv) < 
X^3vol(Si)/r = 3 vo1(S)/t = O(RVr). □ 



Lemma 2.6. Let U be any unit ball whose center is in 
C and at distance 2/3 from 9C. Then U contains 0(1) 
surface area of C. 

Proof: Without loss of generality, assume that U is 
centered at the origin and that (0,0,2/3) is the closest 
point of 9 C to the origin. Let U' be the open ball of 
radius 6 centered at the origin, let V be the open unit 
ball centered at (0, 0, 5/3), and let W be the cone whose 
apex is the origin and whose base is the circle 9U n 9V. 
See Figure 5. U' lies entirely inside C, and since r > 1 , 
we easily observe that V lies entirely outside C. Thus, 
the surface area of 9CnW C 9CnU is at least the area of 
the spherical cap 9U' n W, which is exactly 47t/27. □ 




Figure 5. Proof of Lemma 2.6 



Theorem 2.7. Let S be a set of points in IR^ whose 
closest pair is at distance 1, and let r be any real 
number. Any point in S has O(r^) Delaunay neighbors 
at distance at most r. 

Proof: Let o be an arbitrary point in S, and let B 
be a ball of radius r centered at o. Call a Delaunay 
neighbor of o a friend if it lies inside B, and call a 
friend q interesting if there is another point p G S (not 
necessarily a Delaunay neighbor of o) such that |op| < 
|oq| and z^poq < 1/r. A simple packing argument shows 
that o has at most O(r^) boring friends. 

Let Q be the set of interesting friends of o. Every 
point q G Q lies on the boundary of a Delaunay ball dq 
that contains no points of S in its interior and also has 
o on its boundary. It is straightforward to prove that 
because q is interesting and has distance at least 2 from 
any other point, dq must have radius at least r. Let bq 
be the ball concentric with dq with radius 2/3 less than 
the radius of dq . Finally, for any point q, let Uq be the 
unit-radius ball centered at q . 

We now have a set of unit balls, one for each 
interesting friend of o, whose centers lie at distance 
exactly 2/3 from the boundary of the Swiss cheese C = 
B \ UqeQ • Lemma 2.5, C has surface area O(r^), 
and by Lemma 2.6, each unit ball Uq contains £1(1) 
surface area of C. Since the unit balls are disjoint, it 
follows that o has at most O(r^) interesting friends. □ 
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Theorem 2.8. Let S be a set of points in IR^ whose 
closest pair is at distance 2 and whose diameter is 2A, 
and Jet r he any reed number. There are 0[A^/r] points 
in S with a Delaunay neighbor at distance at least r. 

Proof: Call a point far-reaching if it has a Delaunay 
neighbor at distance at least r, and let Q be the set 
of far-reaching points. Let B be a ball of radius 2A 
containing S. For each q G Q, let fq be a maximal 
empty ball containing q and its furthest Delaunay 
neighbor, and let bq be the concentric ball with radius 
2/3 smaller than fp. By construction, each ball bq has 
radius at least t/2 — 2/3. Finally, for any far-reaching 
point q, let Uq be the unit-radius ball centered at q. 
By Lemma 2.5, the Swiss cheese C = B \ UqeQ ^1 
surface area O(AVr), and by Lemma 2.6, each unit 
ball Uq contains £i(1) surface area of C. Since these 
unit balls are disjoint, there are at most 0(A^/r) of 
them. □ 

Corollary 2.9. Let S be a set of points in IR^ with 
spread A. The Delaunay triangulation of S has com- 
pJexitj 0(A4). 

Proof: For all r, let F(t) be the number of far-reaching 
points in S, i.e., those with Delaunay edges of length 
at least r. From Theorem 2.8, we have F(r) — 0(A^/r). 
By Theorem 2.7, if the farthest neighbor of a point p 
is at distance between r and r + 1, then p has O(r^) 
neighbors. Thus, the total number of Delaunay edges is 
at most 

A A 

^ 0(r2) . (F(r) - F(r + D) = 0(r) • F(r) 

r=0 r=0 

A 

= £o(a3) 

r=0 

= 0(A4) □ 

2.3 Conjectured Upper Bounds 

I conjecture that the lower bounds in Theorem 2.4 
are tight, but Corollary 2.9 is the best upper bound 
known. Nearly matching upper bounds could be 
derived from the following conjecture using a divide and 
conquer argument, suggested by Edgar Ramos (personal 
communication) . 

Let S be a well-separated set of points with closest 
pair distance 1 , lying in two balls of radius A that are 
separated by distance at least cA for some constant 
c > 1 . Call an edge in the Delaunay triangulation of S 
a crossing edge if it has one endpoint in each ball. 



Conjecture 2.10. Some point in S is an endpoint of 
O (A) crossing edges. 

Lemma 2.11. Conjecture 2.10 implies that the Delau- 
nay triangulation of S has 0(min{A^, An, n^)) crossing 
edges. 

Proof: Theorem 2.8 implies that only O(A^) points can 
be endpoints of crossing edges. Thus, we can assume 
without loss of generality that n = 0(A'^). 

We compute the total number of crossing edges by 
iteratively removing the point with the fewest crossing 
edges and retriangulating the resulting hole, say by 
incremental flipping. Conjecture 2.10 implies that we 
delete only 0(A) crossing edges with each point, so 
altogether we delete O(nA) = O(A^) crossing edges. 
Not all of these edges Eire in the original Delaunay 
triangulation, but that only helps us. □ 

Theorem 2.12. Conjecture 2.10 implies that the Delau- 
nay triangulation of n points in JR^ with spread A has 
complexity 0(min{A^ log A, nA, n^j) . 

Proof: Assume Conjecture 2.10 is true, and let S be an 
arbitrary set of n points with diameter A, where the 
closest pair of points is at unit distance. S is contained 
in an axis-parallel cube C of width A. We construct 
a well-separated pair decomposition of S [13], based 
on a simple octtree decomposition of C. The octtree 
has O(logA) levels. At each level i, there are 8^ cells, 
each a cube of width A/2^. Our well-separated pair 
decomposition includes, for each level i, the points in 
any pair of level-i cells separated by a distance between 
cA/2^ and 2cA/2^. A simple packing argument implies 
that any cell in the octtree is paired with 0(1) other 
cells, all at the same level, and so any point appears in 
O(logA) subset pairs. Every Delaunay edge of S is a 
crossing edge for some well-separated pair of cells. 

Lemma 2.11 implies that the points in any well- 
separated pair of level-i cells have 0(A^/8^) crossing 
Delaunay edges. Since there are 0(8^) such pairs, the 
total number of crossing edges between level-i cells is 
O(A^). Thus, there are O(A^logA) Delaunay edges 
altogether. 

Lemma 2.11 also implies that for any well-separated 
pair of level-i cells, the average number of crossing edges 
per point is 0(A/2^). Since every point belongs to a 
constant number of subset pairs at each level, the total 
number of crossing edges at level i is 0(nA/2^). Thus, 
the total number of Delaunay edges is O(nA). □ 

This upper bound is still a logarithmic factor away 
from our lower bound construction when A — o[^/n). 
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However, our argument is quite conservative; all cross- 
ing edges for a well-separated pair of subsets are 
counted, even though some or all of these edges may be 
blocked by other points in S. A more careful analysis 
would probably eliminate the final logarithmic factor. 

3 Nice Surface Data 

Let L be a smooth surface without boundary in IR^. 
The medial axis of L is the closure of the set of points 
in IR^ that have more than one nearest neighbor on L. 
The local feature size of a point x G Z, denoted Ifs(x), 
is the distance from x to the medial axis of L. Let S 
be a set of sample points on I. Following Amenta 
and Bern [1], we say that S is an i-sample of L if the 
distance from any point x G 1 to the nearest sample 
point is at most e • Ifs(x). 

The first step in several surface reconstruction algo- 
rithms is to construct the Delaunay triangulation or 
Voronoi diagram of the sample points. Edelsbrunner 
and Miicke [22] and Bajaj et al. [7, 9] describe algo- 
rithms based on alpha shapes, which are subcomplexes 
of the Delaunay triangulation; see also [25]. Extending 
earlier work on planar curve reconstruction [2, 23], 
Amenta and Bern [1, 3] developed an algorithm to 
extract a certain manifold subcomplex of the Delaunay 
triangulation, called the crust. Amenta et al. [4] 
simplified the crust algorithm and proved that if S is 
an £-sample of a smooth surface L, for some sufficiently 
small £, then the crust is homeomorphic to L. Bois- 
sonnat and Cazals [12] and Hiyoshi and Sugihara [26] 
proposed algorithms to produce a smooth surface using 
natural coordinates, which are defined and computed 
using the Voronoi diagram of the sample points. Further 
examples can be found in [5, 6, 11, 15]. 

In this section, we show that e-samples of smooth 
surfaces can have complicated Delaunay triangulations, 
implying that all these surface reconstruction algo- 
rithms can take quadratic time in the worst case. We 
will analyze our constructions in terms of the sample 
measure of a surface Z, which we define as follows: 



Lemma 3.1. For all e < 1/2, every c-sample of L 
contains D.(]ji{L)/t^] points. 

Proof: Let S be an arbitrary e-sample of Z. Amenta 
and Bern [1] observed that |lfs(p) — lfs(q)| < |pq| 
for any points p,q G L. This observation implies 
that for any point x G I, we have |xp| < £lfs(x) < 



£(lfs(p) + Ixpl] < lfs{p), where p G S is the sample 
point closest to x. Thus, we can cover L with circular 
neighborhoods of radius Ifs(p) around each sample 
point p G S. By similar arguments, the neighborhood 
of p has area at least n{ j^[] — lfs{p))^, and any 
point in the neighborhood of p has local feature size at 
most (1 + ) Ifs(p). It follows that each neighborhood 
has sample measure D.(e^), and since there are n such 
neighborhoods, \i[L) — 0(n£^). □ 

We say that an e-sample is parsimonious if it 
contains 0(|j.(L)/£^) points. 

3.1 Oversampling Is Bad 

The easiest method to produce a surface sample with 
high Delaunay complexity is oversampling, where some 
region of the surface contains many more points than 
necessary. In fact, the only surface where oversam- 
pling cannot produce a quadratic-complexity Delaunay 
triangulation is the sphere, even if we only consider 
parsimonious samples. 

Theorem 3.2. For any smooth non-spherical surface L, 
any £ > 0, and any sufficiently large n, there is a 
parsimonious c-sample of L of size n whose Delaunay 
triangulation has complexity n(rL^). 

Proof: Let S be any parsimonious £-sample of L. Let ct 
be a small sphere intersecting Z in a non-planar curve, 
where the distance from cr to any point os S is at elast 
the radius of cr. Such a sphere always exists unless L is 
itself a sphere. Let a and |3 be extremely short segments 
of the intersection curve I n ff that approximate skew 
line segments. Straighten these curves slightly, keeping 
them on the surface L and keeping the endpoints fixed, 
to obtain curves a' and |3'. Finally, let A and B be sets 
of |S| evenly spaced points on a' and |3', respectively. 
See Figure 6. 




Figure 6. Parsimoniously oversampling a non-spherical surface. 

The Delaunay triangulation of A U B has complexity 
n(|Sp); every point in A is a Delaunay neighbor of 
every point in B. Moreover, any Delaunay circumsphere 
of A U B closely approximates the sphere g and thus 
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excludes every point in S. Thus, S U A U B is a 
parsimonious £-sample of L consisting of n — 3|S 
points whose Delaunay triangulation has complexity 

D.(n^). □ 

The reconstruction algorithm of Amenta et al. [4] 
extracts a surface from a subset of the Delaunay tri- 
angles of the sample points. Their algorithm estimates 
the surface normal at each smple point p using the 
Voronoi diagram of the samples. The cocone at p is 
the complement of a very wide double cone whose apex 
is p and whose axis is the estimated normal vector at p. 
The algorithm extracts the Delaunay triangles whose 
dual Delaunay edges intersect the cocones of all three 
of its vertices, and then extracts a manifold surface 
from those cocone triangles. Usually only a small subset 
of the Delaunay triangles pass this filtering phase, but 
our construction shows that there can be n(n^) cocone 
triangles in the worst case. 

3.2 Uniform Sampling Can Still Be Bad 

Unfortunately, oversampling is not the only way to get 
quadratic Delaunay triangulations. Let S be a set of 
sample points on the surface L. We define the second 
sampling density of a point x G L, denoted sdalx), 
as the distance from x to the second closest sample 
point, divided by Ifs(x). We say that S is a uniform 
t-sample of I if e/4 < sdilx) < e for all x € L} 
Uniform e-samples are also parsimonious e-samples, but 
with absolutely no oversampling. In particular, the size 
of any uniform e-sample is Q[[i(L)/t^). 

Lemma 3.3. For any n and £ > \/T/ri, there is a two- 
companent surface L and an n-point uniform ^-sample S 
of L, such that the Delaunay triangulation of S has 
complexity Cl(n^t^). 

Proof: The surface L is the boundary of two sausages 
Lx and Zy, each of which is the Minkowski sum of a 
unit sphere and a line segment. Specifically, let 

Ix = U + (-W, 0, d + 1 ), (w, 0, d + 1 ) and 

= U+(0,-w,-d-l),(0,w,-d-l), 

where U is the unit ball centered at the origin, w = ne^, 
and d — 4w/£ — 4ne. The local feature size of every 
point on Z is 1 , so any uniform e-sample of Z has 

G((w + ] ]/t^) = e(n) points. 

'There is nothing special about the number 4 here; any 
constant c > 2 will do. However, as c approaches 2, the maximum 
£ for which a c-uniform e-sample exists approaches zero. 



Define the seams and Gy as the maximal line 
segments in each sausage closest to the xy-plane: 

cTx = (— w, 0, d), (w, 0, d) and 
Gy = (0,-w,-d), {0,w,-d). 

Our uniform £-sample S contains 2w/ e -|- 1 points along 
each seam: 

Pi — (i£,0, d) for all integers — w/e < i < w/e, and 
qj — (O.je,— d) for all integers —w/e < j < w/e. 

The Delaunay triangulation of these 0(w/e] — 0(rLe) 
points has complexity 0(w^/e^) = O(n^e^). 

Let Yij be the ball whose boundary passes through 
Pi and qj and is tangent to both seams. This ball may 
contain other portions of the surface, but we claim that 
the intersection is small enough that we can avoid it 
with our sample points. The intersection of Zx and yij 
is a small oval, tangent to Pi and symmetric about the 
plane x = ie. Tedious calculation (which we omit) 
implies that the width of the oval is 

_i / 4dje \ 4w _ 

See Figure 7. 




Figure 7. Two sausages and a sphere tangent to both seams. 

So Zx n Yij lies entirely within a strip of width 2e 
centered along the seam cTx- A symmetric argument 
gives the analogous result for Zy Hyij . We can uniformly 
sample Z so that no other sample point lies within either 
strip. Each segment pTq]" is an edge in the Delaunay 
triangulation of the sample, and there are n(w^/£^) — 
D{n^e^) such segments. □ 

Theorem 3.4. For any n and any £ > ^( log n)/n, 
there is a connected surface Z and an n-point uniform 
e-sample S of Z, such that the Delaunay triangulation 
of S has complexity n(n^£^). 
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Proof: Intuitively, we produce the surface L by pushing 
two sausages into a spherical balloon. These sausages 
create a pair of conical wedges inside the balloon whose 
seams lie along two skew lines. The local feature size 
is small near the seams and drops off quickly elsewhere, 
so a large fraction of the points in any uniform sample 
must lie near the seams. We construct a particular 
sample with points exactly along the seams that form 
a quadratic-complexity triangulation, similarly to our 
earlier sausage construction. Our construction relies 
on several parameters: the radius R of the spherical 
balloon, the width w and height h. of the wedges, and 
the distance d between the seams. 




Figure 8. A smooth surface with a bad uniform £-sample, and a 
closeup of one of its wedges. 

Each wedge is the Minkowski sum of a unit sphere, 
a right circular cone with height h. centered along the 
z-axis, and a line segment of length w parallel to one 
of the other coordinate axes. The boundary of each 
wedge can be decomposed into cylindrical, spherical, 
conical, and planar facets. The cylindrical and spherical 
facets constitute the blade of the wedge, and the seam 
of the blade is the line segment of length w that bisects 
the cylindrical facet. The local feature size of any 
point on the blade is exactly 1 , and the local feature 
size of any other boundary point is its distance from 
the blade. Straightforward calculations imply that the 
sample measure of the wedge is O (w + log h. + 1 ) . 

A first approximation L of the surface I is obtained 
by removing two wedges from a ball of radius R centered 
at the origin. One wedge points into the ball from 
below; its seam is parallel to the x-axis and is centered 
at the point (0, 0, — R + h). The other wedge points into 
the ball from above; its seam is parallel to the y-axis and 
is centered at (0, 0, R — h). Let d = 2R — 2h. — 2 denote 
the distance between the wedges. Our construction has 
1 < w < d < h, so R < 3h.. 

To obtain the final smooth surface L, we round off 
the sharp edges by rolling a ball of radius h./4 inside L 
along the wedge/balloon intersection curves. We call 
the resulting warped toroidal patches the sleeves. The 
local feature size of any point on the sleeves or on the 



balloon is at least h./4. Since L is star-shaped and 
contained in a sphere of radius R, its surface area is 
at most 47tR^ < 367Th.^. It follows that the sleeves 
have constant sample measure. The local feature size of 
wedge points changes only far from the blades and by 
only a small constant factor, so ]i{L) = 0(w+logh.-|-l ). 
To complete the construction, we set w = ne^, d — 4n£, 
and h, = lOne. See Figure 8. 

Finally, we construct a uniform £-sample S with 
©[w/e] sample points evenly spaced along each seam 
and every other point at least £ away from the seams. 
Setting h. > 5d (and thus R > lOd) ensures that the 
Delaunay spheres ytj between seam points do not touch 
the surface except on the blades. By the argument in 
Lemma 3.3, there are Q.(w^/e^) =a{n^£2) Delaunay 
edges between seam points. □ 

3.3 Some Surfaces Are Just Evil 

In this section, we describe a family of surfaces for 
which any parsimonious £-sample has a Delaunay trian- 
gulation of near-quadratic complexity. First we give a 
nearly trivial construction of a bad surface with several 
components, and then we join these components into 
a single connected surface using a method similar to 
Theorem 3.4. 

Lemma 3.5. For any n and any e < -^/l/n, there is a 
smooth surface L such that the Delaunay triangulation 
of any parsimonious e-sample of L has complexity 
nin^E"*), where n is the size of the sample. 

Proof: Let P be a set containing the following k points: 

= (ik, 0, k^) for all integers -k/4 < 1 < k/4, and 
qj = (0, jk, -k^) for all integers -k/4 < j < k/4. 

We easily verify that every pair of points pi and qj lie 
on a sphere ytj with every other point in P at least unit 
distance outside. 

Let L — Upep > where Up is the unit-radius sphere 
centered at p. Clearly, Ifs(x) — 1 for every point x G L, 
so |j.(L) — 47Tk. Let S be an arbitrary parsimonious e- 
sample of Z, let n = jS| = 0(k/£^), and for any point 
p e P, let Sp = S n Up be the sample points on its unit 
sphere. 

Choose an arbitrary Delaunay pair Pi,qj € P, and 
let Y be a sphere concentric with ytj but with radius 
smaller by 1. This sphere is tangent to Up^ and Uq^ 
but is at least unit distance from every other component 
of L. Expand y about its center until it hits (without 
loss of generality) a point p' e Sp^ , and then expand it 
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about p' until it hits a point q' G Sp^. The resulting 
sphere y passes through p' and q' and has no points 
of S in its interior, so p' and q' are joined by an edge 
in the Delaunay triangulation of S. There are at least 
a(k2) = Oln^e^) such edges. □ 

To create a connected surface where good sample 
has a complicated Delaunay triangulation, we add 
'teeth' to our earlier balloon and wedge construction. 
Unfortunately, in the process, we lose a logarithmic 
factor in the Delaunay complexity. 



Theorem 3.6. For any n and any e < ^ (log n)/n, 
there is a smooth connected surface L such that the 
Delaunay triangulation of any parsimonious c-sample 
ofL has complexity D.(n^ e'^ / log^ (nc^)) , wheren is the 
size of the sample. 

Proof: Intuitively, we create the surface L by pushing 
two rows of regularly spaced unit balls into a large 
spherical balloon, similarly to the proof of Theorem 3.4. 
As before, the surface contains two wedges, but now 
each wedge has a row of smeJl conical teeth. Our 
construction relies on the same parameters R,w,h of 
our earlier construction. We now have additional 
parameter t, which is simultaneously the height of the 
teeth, the distance between the teeth, and half the 
thickness of the 'blade' of the wedge. 

Our construction starts with the (toothless) surface 
described in the proof of Theorem 3.4, but using a ball 
of radius t instead of a unit ball to define the wedges. 
We add w/t evenly-spaced teeth along the blade of 
each wedge, where each tooth is the Minkowski sum 
of a unit ball with a right circular cone of radius t. 
Each tooth is tangent to both planar facets of its 
wedge. To create the final smooth surface L, we roll 
a ball of radius t/3 over the blade/tooth intersection 
curves. The complete surface has sample measure 
e((w/t)(l + logt) + logh+1). Finally, we set the 
parameters w = t^, h = t^, and R = 20t^, so that 
^(I) =0(tlogt). 

Let S be a parsimonious e-sample of Z, and let u = 
\S\ = 0((tlogt)/e'^). For any pair of teeth, one on each 
wedge, there is a sphere tangent to the ends of the teeth 
that has distance £1(1 ) from the rest of the surface. We 
can expand this sphere so that it passes through one 
point on each tooth and excludes the rest of the points. 
Thus, the Delaunay triangulation of S has complexity 
a(t2) =a(n2e4/log^(n£2)). □ 



3.4 Randomness Doesn't Help Much 

Golin and Na recently proved that if S is a random set 
of n points on the surface of a convex polytope, then the 
expected complexity of the Delaunay triangulation of S 
is 0(n) [24]. Unfortunately, this result does not extend 
to nonconvex objects, even the random distribution of 
the points is proportional to the sample measure. 

Theorem 3.7. For any n, there is a smooth connected 

surface L, such that the Delaunay triangulation of n 
independent uniformly-distributed random points in L 
has complexity Q[n^ / log n] with high probability. 

Proof: Consider the surface I consisting of 0(n/ logn) 
unit balls evenly spaced along two skew line segments, 
exactly as in the proof of Theorem 3.5, with thin cylin- 
ders joining them into a single connected surface. With 
high probability, a random sample of n points contains 
at least one point on each ball, on the side facing the 
opposite segment. Thus, with high probability, there 
is at least one Delaunay edge between any ball on one 
segment and any ball on the other segment. □ 

Theorem 3.8. For any n, there is a smooth connected 
surface L, such that the Delaunay triangulation of 
n independent random points in L, distributed pro- 
portionally to the sample measure, has complexity 
6(n^/ log'* n] with high probability. 

Proof: Let L be the surface used to prove Theorem 3.6, 

but with 0(u/log tl] teeth. With high probability, a 
weighted random sample of L contains at least one point 
at the tip of each tooth. □ 

4 Conclusions 

We have derived new upper and lower bounds on the 
complexity of Delaunay triangulations under two dif- 
ferent geometric constraints: point sets with sublinear 
spread and good samples of smooth surfaces. Our 
results imply that with very strong restrictions on the 
inputs, existing surface reconstruction algorithms are 
inefficient in the worst case. 

Our results suggest several open problems, the most 
obvious of which is to tighten the spread-based bounds. 
Even the special case A = 0(n^/^) is open. 

Another natural open problem is to generalize 
our analysis to higher dimensions. Using the proof 
techniques in Section 2.2, we can show that any 
d-dimensional Delaunay triangulation has 0(A'^+^) 
edges. We conjecture that the total complexity is 
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always 0(A'^) and can only reach the maximum 

n(nrd/21) when A = Q(u). 

Our bad surface examples are admittedly contrived, 
since they have areas of very high curvature relative 

to their diameter. An interesting open problem is 
whether there are bad surfaces with smaller 'spread', 
i.e., ratio between diameter and minimum local feature 
size. What is the worst-case complexity of the Delaunay 
triangulation of good surface as a function of the spread 
and sample measure of the surface? 

Our results imply that any Delaunay-based surface 
reconstruction algorithm can be forced to take super- 
linear time, even for very natural surface data. It may 
be possible to improve these algorithms by adding a 
small number of Steiner points in a preprocessing phase 
to reduce the complexity of the Delaunay triangulation. 
In most of our bad surface examples, a single Steiner 
point reduces the Delaunay complexity to 0(n). Bern, 
Eppstein, and Gilbert [10] show that any Delaunay 
triangulation can be reduced to 0(n) complexity in 
O(nlogn) time by adding O(rL) Steiner points; see 
also [14]. Unfortunately, the Steiner points they choose 
(the vertices of an octtree) may make reconstruction 
impossible. In order to be usable, any new Steiner 
points must either lie very close to or very far from the 
surface, and as our bad examples demonstrate, both 
types of Steiner points may be necessary. Boisson- 
nat and Cazals (personal communication) report that 
adding a small subset of the original Voronoi vertices as 
Steiner points can significantly reduce the complexity 
of the resulting Voronoi diagram with only minimal 
changes to the smooth surface constructed by their 
algorithm [12]. 

Very recently, Dey et al. [18] developed a surface 
reconstruction algorithm that does not construct the 
entire Delaunay triangulation. Their algorithm runs in 
O(nlogn) time if (loosely speaking) the density of the 
sample points varies smoothly over the surface. 

Finally, are there other natural geometric conditions 
under which the Delaunay triangulation provably has 
small complexity? 
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